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Abstract 



We use path-integrals to derive a general expression for the semiclassi- 
cal approximation to the partition function of a one-dimensional quantum- 
mechanical system. Our expression depends solely on ordinary integrals which 
involve the potential. For high temperatures, the semiclassical expression is 
dominated by single closed paths. As we lower the temperature, new closed 
paths may appear, including tunneling paths. The transition from single to 
multiple-path regime corresponds to well-defined catastrophes. Tunneling sets 
in whenever they occur. Our formula fully accounts for this feature. 
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We may use path-integrals in the description of Quantum Statistical Mechanics. The 
partition function for a one- dimensional quantum-mechanical system, for instance, can be 
expressed as 

poo px{/3h)=x 

Z((3) = dx [Dx(t)} e~ s/h , (la) 

J — oo J x(0)=xq 



ppn 
S= dr 



o 



\m± 2 + V(x) 



(lb) 



For an arbitrary potential, V(x), this integral has been approximated by means of pertur- 
bative variational |l|J^] and numerical techniques. Here, we shall concentrate on 

the semiclassical approximation to the integral in order to: (i) derive a general formula for 
Z sc (f5), for an arbitrary potential, that does not require a detailed knowledge of the classical 
motion;[] (ii) discuss the onset of tunneling and relate it to the study of singularities in 
Catastrophe Theory. 

The semiclassical evaluation of (|la|) yields: 

Z sc ((3)= dx £ e- s ^^ h A- 1/2 (x ,f3), (2) 

J-oc j=1 

where Sj(xo,/3) denotes the action and Aj(xo,P) represents the determinant of the fluctua- 
tion operator, 

F = -M^ + V[x j (T)], (3) 

both calculated at the j-th classical trajectory, Xj(r), satisfying the boundary conditions 
Xj(0) = Xj(/3h) = xo] N(xo,f3) is the number of classical trajectories which are minima of 
the action functional. (Note that Xj(r) are closed paths, but not necessarily periodic. If the 
integral in (0) is evaluated using the saddle-point approximation, then only the latter must 
be considered [|J. However, there are situations in which this further approximation cannot 
be used, e.g., if V(x) = Ax 4 .) 



1 By classical motion we mean motion satisfying the Euler-Lagrange equation Mx — V'(x) 
which is the equation of motion for a particle moving in a potential minus V. 
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The action has a simple expression in terms of the turning points of the classical trajec- 
tory: 

3 3 
P X i . P X 

Sj(x ,P) = (3hV(x>±) ±2 dxMv(x,x j ±) + 2n / . dxMv{x,xi) . (4) 

Jxq J x 3 _ 



Here, v(x,y) = y (2/M)[V(x) — V(y)} and x 1 + (x 3 _) are turning points to the right (left) 
of xq. The first term in corresponds to the high-temperature limit of Z(/3), where the 
classical paths collapse to a point, i.e., x± — > xq. For motion in regions where the inverted 
potential is unbounded (hereafter called unbounded motion), n — 0. For periodic motion, n 
counts the number of periods and the second term is absent. For bounded aperiodic motion, 
all terms are present. The last two terms will be negligible for potentials which vary little 



over a thermal wavelength, A = KyP/M. However, by decreasing the temperature, quantum 
effects become important. 

For trajectories having a single turning point (n = 0), x± are given implicitly by 

f x ± dx , , 

Ph = ±2 ~i jt j (5) 

Jx V{X,X±) 



and the fluctuation determinant by 



A 1 MV'(x 3 ± ) 



8(3 



dxi 



± 



(6) 



The formula above is only valid for trajectories with a single turning point (n = 0). For 
those with two turning points, we have derived a generalization which depends on the 
number of periods, n. However, the latter are naturally excluded from the calculation, as 
we shall see in the sequel. 

Formulae (|) and @ do not require full knowledge of the classical trajectory, as the 
dependence on Xj(r) comes only through the turning point. To see how this can simplify 
the evaluation of Z sc (/3), let us take the harmonic oscillator, V(x) = \Moj 2 X . cLS clll example. 
In this case, given xq and f3 there is only one trajectory, with a single turning point, given 
by x + (x_) = Xq/ cosh(/3hui/2) for x < (> 0). S and A can also be readily calculated, the 
final result being 



oo 



At „-{Muxl/h)tank{Phu/2) 



Mu 



(7) 



— oo 



2tt% smh{f3huj) ' 



which in this case is exact. 



For V(x) = ^Mu 2 x 2 + Ax 4 , A > 0, it is still true that there is only one trajectory, given 
x and P, with a single turning point. Eq. (|5]) yields: 



x± 



x cn.(u±,k±), 



(8a) 



(3hu / AXx 2 ± 




(8b) 



Muj 2 + 2Xx 2 ± 
Mu 2 + 4\x 2 ± ' 



(8c) 



where cn(u, k) is one of the Jacobian elliptic functions ||. We may, then, use Eqs. (f|) and 
and change the integration variable in (0) from xq to x± to obtain Z sc {(3). The resulting 
integral can be evaluated numerically. 

A more interesting situation occurs in the case of the anharmonic oscillator, V(x) = 
X(x 2 — a 2 ) 2 . For x 2 > a 2 only single paths with single turning points exist for fixed Xq and 
(3. However, there is also a region, x 2 < a 2 , where the classical motion is bounded and a 
much richer structure exists |7j], one in which more than one classical path may exist for 
given values of xq and (3. 

In a region of bounded classical motion (a well in — V), such as x 2 < a 2 for the anharmonic 
oscillator, the number of classical trajectories changes as the temperature drops. If < /3 < 



7t/hu) m (where u m = y— V"{x m )/M and x m is a local minimum of — V), for every xq in 
this region there is only one closed path, with a single turning point, satisfying the classical 
equations of motion. For x < x m (> x m ) this path goes to the left (right) and returns to 
xq. For xq = x m , it sits still at the bottom of the well. It is this single-path regime which 
goes smoothly into the high-temperature limit. 

For (3 = ii/huj m , the solution that sits still at x m becomes unstable. Its fluctuation 
operator is that of a harmonic oscillator with u 2 = uj 2 m . Its finite temperature determinant 
is A(x m ,/3) = 2nhsva.{j3huj m ) / Mu m . This goes through zero at (3 — 7i/hu m and becomes 
negative for (3 > n/hum, thus signaling that x{r) = x m becomes unstable. At the same 
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time, two new classical paths appear, symmetric with respect to x m , as depicted in Fig. 1. 
Therefore, at Xq = x m , we go from a single-path regime to a triple-path regime as we cross 
(3 = n/huj m . The two new paths are degenerate minima, whereas the path that sits still at 
x m becomes a saddle-point of the action, with a single negative mode. 

As (3 grows beyond 7r/hu m , an analogous situation occurs for other values of x inside the 
well. When the fluctuation determinant around the classical path for a given xo vanishes, a 
new classical path appears. Its single turning point lies opposite, with respect to x m , to that 
of the formerly unique path. It may be interpreted as a tunneling trajectory, since it traverses 
a classically forbidden region of V. If (3 is increased further, this tunneling trajectory splits 
into two, as illustrated in Fig. 2. One is a local minimum of the action, whereas the other is 
a saddle-point, with only one negative mode. Again, we have transitioned from a single to a 
triple-path regime. As (3 grows, the triple-path region spreads out around x m . The frontiers 
of that region are defined by the points xo such that (8(3 /dx±) = 0. 

The phenomenon just described is an example of catastrophe. It takes place whenever 
the number of classical trajectories changes, with one or more of them becoming unstable, as 
we lower the temperature. Conversely, we may say that the phenomenon is characterized by 
the coalescence of two or more of the classical trajectories, as we increase the temperature. 
This is akin to the ocurrence of caustics in Optics [§], where light rays play the role of 
classical trajectories and the action is replaced with the optical distance. 

If we denote by S\ the coordinate associated with the direction of instability of our exam- 
ple, the action can be viewed as S = S(si, . . . ; x , (3). If we use the basis of eigenf unctions 
of the fluctuation operator around classical paths, S\ will correspond to the direction along 
the one eigenfunction whose eigenvalue goes through zero, with the dots referring to all 



others. Catastrophe Theory p|-|T0f allows us to write down the "normal form", Sn, of this 
generating function in the three-dimensional subspace made up by the unstable direction 
of state space (the set of paths) and the two control variables (related to (3 and xq); it is 
this subspace which is relevant for the study of the onset of instability. As catastrophes are 
classified by their codimension in control space, which here is two-dimensional, only those of 



codimension 1 (the fold) or 2 (the cusp) are generic (i.e., structurally stable). The pattern 
of extrema then leads to the cusp, whose generating function is 

SW(si;ui,Vi) = sf + msl + VxS! , (9) 

where U\ and Vi, the control parameters, are related to j3 and Xq, respectively. This catas- 
trophe is defined by 



dS N d 2 S N 



0. (10) 



dsi ds\ 

Eliminating s\ from these equations, we can draw the bifurcation set in control space (the 
curve v\ = — 4uf/27), as well as the pattern of extrema of SV (Fig. 3). Assuming that u\ and 
v\ are linearly related to {J5 — n/hu m ) and (xq — x m ), respectively, for small values of both, 
we can also schematically plot the classical action (i.e., the action for classical trajectories) 
as a function of xq for different values of f3%\ see Fig. 4(a,b). 

New tunnelings, accompanied by new catastrophes, will occur as we keep increasing f3. 
From (3 = 2^/huj m , we start having tunneling trajectories with x(0) = 0, i.e., with one 
full period. For these, the determinant of fluctuations vanishes because the first, not the 
second bracket of (H) goes trough zero. However, the fluctuation operator around such 
tunneling trajectories already has a negative eigenmode. This follows from the fact that 
the zero-mode, given by x c \(t), has a node. Thus, it cannot be the ground-state for the 
associated Schrodinger problem. This means that a new catastrophe takes place along a 
second direction in function space. 



2 This plot can be obtained by exploiting the relation between the cusp and the swallow- 
tail catastrophes. The latter has a generating function whose normal form is W(si;a,b,c) = 
sf + asf + bs\ + csi, where a, b, c are control variables. The extremum condition is, then, 
dW/dsi = 0. The identifications £V = — c/5, u\ = 3a/5 and v\ = 26/5 will make the additional 
condition that defines the swallowtail, d 2 W/ds\ = 0, coincide with the requirement OSn/Qsi = 0. 
(In the jargon of catastrophe theory: the bifurcation set for the swallowtail coincides with the 
equilibrium hyper surf ace of the cusp.) 
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This new catastrophe is also a cusp. At x m , as (3 becomes larger than 2 r K/Tmo m) the 
solution that sits still acquires a second negative eigenvalue - - thus becoming unstable 
along a new direction in function space — and two other solutions appear. They both 
have a negative eigenvalue along the first direction of instability and a positive eigenvalue 
along the new direction of instability. If we combine the information from both directions, 
we find that we go from two minima and a saddle-point with only one negative eigenvalue 
(hereafter, a one-saddle) to two minima, two one-saddles and one two-saddle. The two one- 
saddles are degenerate in action and represent time-reversed periodic paths. As (3 increases 
beyond 27r/hu> m , the same phenomenon takes place for xq around x m : the one-saddles that 
already existed for (3 < 2n/huj m in the three-path region around x m become unstable along 
a second direction, and a five-path region grows around x m , with two minima, two (periodic, 
degenerate in action) one-saddles and a two-saddle. 

The new cusp can be cast into normal form using the second direction of instability: 

S N (s 2 ; u 2 , v 2 = 0) = 4 + u 2 s 2 2 ~ (si + y ) 2 . (11) 

The absence of the linear term in Eq. flTT| ) [compare with Eq. (^)] reflects the degeneracy 
in action of the two one-saddles. However, there is another degeneracy: since V does 
not depend on (3, if x c \(r) is a solution of the Euler-Lagrange equations, so is x c \(r + r ). 
If x c i(t) is periodic, x c \{t + r ) describes the same path (and so has the same action), 
but with another starting point. This can be represented in Eq. (|TTD by the choice u 2 oc 
[(^o - x m ) 2 + k(2ir/hu m - (3)}, with k > 0. See Fig. 4(c). 

As we approach (3 = ?>ir/hLo m) the situation becomes similar to the one near (3 = n/huUm. 
The difference is that we now have to deal with a tunneling path with more than one and less 
than two periods. Near (3 = 4n/tuj m , two-period tunnelings intervene, a situation similar to 
that at (3 = 2ir /Tiuj m . The pattern which develops ]11[ is depicted in Fig. 5. 



Despite the fact that we keep adding new extrema as we lower the temperature, many 
of which represent tunneling solutions, only two of these are stable (i.e., minima). In a 
semiclassical approximation in euclidean time, appropriate to equilibrium situations, these 



are the only extrema we have to sum over, meaning that N(xq, 0) never exceeds two. The 
local minimum is a tunneling solution, whereas the global one becomes the unique solution 
in the high-temperature regime. For (3 < Tr/huj m , there will be regions with either one or 
two minima of the classical action. The transition from a single-minimum to a double- 
minimum region occurs at values of Xq where the fluctuation determinant vanishes due to 
the appearance of a caustic, thus leading to a singularity in the integrand of Eq. (0). This 
is not a disaster, however, as this singularity is integrable. (Such a singularity is an artifact 



of the semiclassical approximation. It disappears in a more refined approximation fq,|12Hl4 



in which one includes higher order fluctuations in the direction(s) of function space where 
the instability sets in.) We shall exhibit the details of this calculation in a forthcoming 
publication. 

To conclude, two remarks: (i) our results may be used to derive a semiclassically improved 
perturbation expansion, which will be presented elsewhere; (ii) the semiclassical partition 
function incorporates, in an approximate way, quantum effects that become more and more 
relevant with the increase of the thermal wavelength. At very low temperatures, however, 
the quadratic approximations inherent to the semiclassical approximation fail to capture 
the details of the potential, which are important in the regime of large thermal wavelengths. 
Thus, we expect our calculations to be a better approximation at high temperatures. Indeed, 
studying the T = limit of Z sc {(3) for the anharmonic oscillator, one finds no corrections 
to the unperturbed ground state energy. Nevertheless, this can be corrected by using the 
semiclassically improved perturbation theory just mentioned. 
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Figure Captions: 



Figure 1: (a) Classical paths at x m for (3 < n/hum (1) and (3 > ir/hu m (2 and 3); (b) 
Sketch of how the extrema change along the unstable direction in function space. 

Figure 2: (a) Classical paths at x m for (3 < (3 C (1), f3 — f3 c (1 and 2&3) and (3 > (3 C (1,2 
and 3), 7i/hu m < (3 C < 2ir/hu m ; (b) Sketch of how the extrema change along the unstable 
direction in function space. 

Figure 3: Bifurcation set for the cusp; pattern of extrema shown schematically. 

Figure 4: Evolution of the classical action as (3 changes, (a) < f3 < n/ficu m ; (b) n /%uo m < 
(3 < 2ii/hu m ; (c) 2ii/hu m < (3 < ?>n/Tiuj m . 

Figure 5: Partition of control space into p-solution regions (p = 1, 3, 5, 7, . . .). 
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